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Information-theoretical quantities such as statistical distinguishability typically result from opti¬ 
misations over all conceivable observables. Physical theories, however, are not generally considered 
valid for all mathematically allowed measurements. For instance, quantum field theories are not 
meant to be correct or even consistent at arbitrarily small lengthscales. A general way of limiting 
such an optimisation to certain observables is to first coarse-grain the states by a quantum channel. 

We show how to calculate contractive quantum information metrics on coarse-grained equilibrium 
states of free bosonic systems (Gaussian states), in directions generated by arbitrary perturbations 
of the Hamiltonian. As an example, we study the Klein-Gordon field. If the phase-space resolution 
is coarse compared to h, the various metrics become equal and the calculations simplify. In that 
context, we compute the scale dependence of the distinguishability of the quartic interaction. 


The application of tools from quantum information 
theory (QIT) to QFT is not an entirely straightforward 
matter. While standard physical applications of a theory 
only requires that one be able to compute expectation 
values of certain specific observables on specific states, 
QIT often requires optimisations over all possible states 
or observables. Hence, it requires a very detailed un¬ 
derstanding of the operational domain of validity of the 
theory under study. 

The main formalism for QFT comes with several 
unique features which can potentially complicate such 
analysis, such as its lack of a Hilbert space or algebraic 
formulation, or the unavoidable use of divergent asymp¬ 
totic series. Moreover, most interacting QFTs are not 
meant to be valid in the continuum. Instead, they are 
defined only relative to an unphysical ultraviolet regula¬ 
tor which can be thought of as mimicking an underlying 
discrete space. The process of renormalisation consists in 
running the parameters of the theory as function of the 
regulator in such a way that predictions are independent 
of it. But even for the simplest interacting QFTs such as 
quantum electrodynamics, this game fails below a certain 
finite lengthscale. 

Certainly, in an optimisation involving in principle all 
observables, one should avoid those for which the the¬ 
ory’s predictions are not trusted. A simple way to deal 
with this issue is to add a cutoff on all momentum in¬ 
tegrals involved in calculations, which usually amounts 
to “tracing-out” the high momentum modes m- Here, 
we want to examine this question more carefully, by ex¬ 
plicitly considering a physical (experimental) limit on the 
resolutions (spatial or otherwise) of the accessible observ¬ 
ables. 

In addition to a limit on spatial resolutions, we study 
the consequence of a finite resolution on the measurement 
of local field values (quadratures). While the spatial res¬ 
olution makes high momenta effectively unobservable, a 
low field value resolution hides high order polynomials 
in the field operators [J]. We show that for the Klein- 
Gordon field, these resolution parameters have a qual¬ 
itative influence on the asymptotic dependence of dis¬ 
tinguishability metrics on scale (which we define as the 


minimal spatial resolution). 

For this purpose, we introduce techniques which al¬ 
low for the calculations of distinguishability measures be¬ 
tween the thermal equilibrium states (or ground states) 
corresponding to perturbations of any given quadratic 
bosonic Hamiltonian, and coarse-grained by any Gaus¬ 
sian channel (to be defined below). Since we only con¬ 
sider the effect of the perturbation on the state to first 
order in the coupling, this amounts to computing the 
components of Riemannian metrics on the manifold of 
Gaussian states, but also in non Gaussian tangent direc¬ 
tions. These calculations can in principle be extended 
outside of the Gaussian manifold perturbatively. 


The main results are formula for the calculation of 
generating functions for the components of the coarse¬ 
grained metrics, which can then be used to compute the 
quantity defined in Equ. (24) (an example is given in 
Equ. (69)). The generating function is computed in full 
generality for the metric (Equ. (65)) and the classical 
Fisher information metric (Equ. (73 )). Moreover, an ap¬ 
proximate formula (Equ. ( |76[ )) is shown to be correct for 
all metrics in the limit of large quadrature imprecision 
for perturbations of a scalar field, where ATapo is given in 

Equ. (ITTTl). 


In the context of QFT, this allows for the calculation 
of the coarse-grained distinguishability between the per¬ 
turbed ground states or thermal states of free fields. Such 
calculations can be directly useful for the characterisation 
of quantum phase transitions mm, and in a recently 
proposed approach to renormalisation HE]. In addition, 
integrating such quantity along Hamiltonian paths can 
be used to obtain bounds on ground-state entanglement, 
as discussed in the outlook. Such calculations are of in¬ 
terest in condensed matter theory, as well as in relation 
to that AdS/GFT correspondance [7]. 


Our use of the Gaussian state formalism is a somewhat 
unusual approach to quantum field theory, however it has 
several advantages: it allows in principle for an exact 
mathematical formulation of free field theory (states on 
a G*-algebra), and it makes it easy to consider the action 
of a large class of quantum channels. 
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I. QUANTUM INFORMATION METRICS 

We focus on the calculation of quantum information 
metrics which directly generalise the classical Fisher in¬ 
formation metric. These are Riemannian metrics on the 
manifold of mixed states. Infinitesimally, these metrics 
measure statistical distinguishability between states, such 
as in the context of hypothesis testing or parameter esti¬ 
mation. 

In a C*-algebraic framework, the main stage is the 
algebra A generalising the set of linear operators on a 
Hilbert space. A state is defined as a positive function p 
which assigns a probability p{A) to any effect A, i.e., any 
self-adjoint element A G A with spectrum between 0 and 
1. In finite dimension p can be represented by a density 
matrix p such that p{A) = "Tv{pA). 

A Riemannian metric is defined by a scalar product 
defined in the tangent space at every point of a manifold. 
A tangent vector to a state p with density matrix p can be 
defined by a traceless self-adjoint operator X (because, 
for e small enough, p-\-eX must be also a density matrix). 
Hence we write the tangent space at p as 

Tp = {Xt=A,Tr(A)=0}. (1) 

The positivity constraint does not appear here because 
we are interested in the bulk of the manifold. We will 
deal with boundaries (where the density matrix is not 
full rank) by approaching them from the bulk. 

A metric associates to every p a positive linear operator 
on Tp defining the scalar product 

^ T^{XVL-p\Y)). ( 2 ) 

The quantity eTr(AQ“^(Ar))5 is the distance between p 
and p + eX as e —?► 0. 

Of particular interest are the Riemannian metrics 
which are contractive under the action of any channel 
(see below). This contractivity is required for the metric 
to represent any type of information-theoretic quantity 
such as distinguishability. In fact, classically, the con¬ 
tractivity condition selects a single metric: the Fisher 
information metric. 

In quantum theory, Petz and Sudar [5] showed that 
these contractive metrics are one-to-one with operator 
monotone function 0 : K+ —>■ K+ such that 9{f) = t9(t~^) 
for all t > 0. An operator monotone function has the 
property that, when applied to operators via functional 
calculus, 9{A) < 9{B) whenever A < B (i.e., B — A is 
positive). 

The function 9 defines the kernel via its inverse 
Up as follows: 

Up = 9{LpR-p^)Rp, (3) 

where Rp{A) := Ap and Lp{A) := pA for any matrix A. 

For instance, in classical probability theory, which is 
equivalent to considering only density matrices diagonal 


in the same basis (the diagonal elements then are a proba¬ 
bility distribution), since superoperators Rp and Lp com¬ 
mute, we obtain simply Op = 9(l)Rp. The number 0(1) 
is an unimportant overall factor which we can pick to be 
0(1) = 1, leading to 

Op (A) = Ap (classical), (4) 

which is independent of the function 0. The resulting 
metric on probability distributions is simply the Fisher 
information metric. Therefore, the contractive metric 
parametrised by 0 are all possible non-commutative gen¬ 
eralisations of the Fisher information metric. 

We will see that it is generally more convenient to work 
in the cotangent space at p, namely the linear dual T* of 
Tp. It is the set of linear functionals on the real vector 
space Tp. These functionals can be characterised by ma¬ 
trices = A through X Tr(AAi). Moreover, since 
Tr(A') = 0, we are free to choose an additional constraint 
which we take to be Tr(A/5) = 0. In this manner, we 
simply have Tr((p -|- eX)A) = Tr(A'A), which directly 
gives an interpretation of the functional associated with 
A as an observable. Moreover, this extra condition makes 
sense in infinite dimension where the trace may not ex¬ 
ists. Hence we write the cotangent space as 

t; = {At = A, p(A) = 0}. (5) 

The metric induces also a scalar product on Tf which 
is {A,B)p = Tr(AQp(i?)). In infinite dimension, we can 
remove the trace by defining 

Qp:=9{LpRf^), (6) 

so that RpQp = Up, leading to the expression 

(A,H)p = p(A0p(R)), (7) 

where we used the fact that [Rp, Lp] = 0. Observe more¬ 
over, that if p is a thermal state p~^^, then the superop¬ 
erator LpRf^ simply yields the imaginary time evolution: 

LpRf^{A) = pAp~^ = e~^^ Ae^^. ( 8 ) 

One can freely move between the tangent 

(Schrodinger) picture and the cotangent (Heisen¬ 

berg) picture by using the metric kernel 0“^. Indeed, 
contracting the metric with a tangent vector gives a 
cotangent vector. Hence, if Qp(A) = X and Ltp{B) = Y, 
then 

Tr(AQ;i(r)) = (A,H)p. (9) 

Particular contractive metrics which appear in the lit¬ 
erature are given by differentiating certain measure of 
distinguishability which have operational interpretations. 
That is, the geodesic distance matches the corresponding 
measure to lowest order (in the distance). 

For instance, differentiating the relative entropy yields 
the Kubo-Mori metrics defined in Equ. (|I^ below. It is 
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of particular importance to us because of its relation to 
first-order perturbation theory, as explained in the next 
section. 

Perhaps the most important metric is the Bures met¬ 
ric, defined by 9{t) = (1 -I- 1)/2 because it has several 
nice features. Its geodesic distance has a closed analyt¬ 
ical form as the Bures distance [5]. It also gives a tight 
bound on the variance of parameter estimation [lOL lllj , 
and is as such usually called the quantum Fisher infor¬ 
mation. Moreover, the Bures metric is the smallest of 
the contractive metrics (normalised by 6*(1) = 1) [5]. 

1 

We will also refer to the metric defined by 9(t) = t‘^, 
because it is especially easy to compute, and is involved 
in a quantum generalisation of the classical test [T^. 

We have 

d{8,{p^,),84p^))=X\\V2-VP\e,p + OiX^) (15) 

with 

\\V\\l^ := /32Tr(£:*H^(H)H-i^£:*H^(y)) (16) 

and Hp(H) is the first-order term in the imaginary time 
Dyson series: 

:= [" p-^Vp^-^ds 

Jo 

(17) 

II. COARSE-GRAINED DISTINGUISHABILITY 

__ 1 d o(^H+\V)\ 

~ pZdX 

In what follows, given a linear “superoperator” 8 de¬ 
fined on the algebra of observables, we write £* for 
its “pre-dual”, namely its adjoint with respect to the 
Hilbert-Schmidt inner product: 

because 81^ (p) = 1) Incidentally, the inverse of the su¬ 
peroperator also dehnes one of the monotone metrics: 

the Kubo-Mori metric, defined from the function 

9^{x) = f x^ds. (18) 

p{£{A)) = Tr{£{A)p) = Tr(A£*(p)) (10) 

Jo 

for any observable A and any state p. We are particu¬ 
larly interested in the case where f is a unital completely 
positive map, that is, £[ 1 ) = 1, and {£ (8)idn)(>l) > 0 for 
all A > 0 and all finite extra dimension n. These condi¬ 
tions guarantee that p o {8 ® idn) (with density matrix 
(5* 0 id)(/5)) is a valid state whenever p is. 

If 8 is such a unital completely positive map, it rep¬ 
resents the action of a quantum channel in the Heisen¬ 
berg picture, while £* represents the same transforma¬ 
tion in the Schrodinger picture. On infinite-dimensional 
systems, 8 is always defined, but may not be. 

Let 

Equ. (16) needs to be reformulated in a way which is 
tractable and suitable for the Gaussian and QFT formal¬ 
ism. As we have seen in the previous section, the metrics 
are easier to formulate in the Heisenberg picture (i.e. on 
the cotangent space). 

For this purpose it proves useful to define the map 

Tip := , (19) 

which is the linearisation of the diffeomorphism H i—>■ 
H' = log5*(e“^^), where H' can be thought of as 
a coarse-grained effective Hamiltonian since = 

8t,{e~d^). Note that the map H H' is the renor- 

(11) 

malisation group transformation introduced in Ref. |13j . 
This allows us to rewrite the norm in Equ. (161 as 

be the thermal state for the Hamiltonian H at inverse 
temperature /3. Geometrically, we want to compute the 
information metrics pulled back to the manifold of Hamil¬ 
tonians. 

Physically, this has the following interpretation. Given 
the metric which gives rise to the geodesic distance 

d, and given the channel 8, and the Hamiltonian H, we 
want to compute the coarse-grained distance 

||Ei,, = /32Tr(£.H^(H)7^p(E)) 

= P^TriVnpTZpiV)) 

where we used the cyclicity of the trace and the fact that 
the metrics are symmetric. 

Writing 

(A,H)p = Tr(AHp(i?)) (21) 

d{84p^,),8M)) (12) 

and 

to lowest order in A where 

(A,H)^ = Tr(AH^(H)) (22) 

(13) 

for the cotangent scalar products associated to the two 
metrics, one can check that TZp introduced above is de- 

are also normalised thermal states. For notational con¬ 
venience, we assume that 

fined by the fact that it is the adjoint of 8 with respect 
to these two scalar products as follows 

Tr(/5(y2 - Vi)) = 0. (14) 

(A,£-(H))f = (7^p(A),H)po£ (23) 


MP{V2 - Vi)) = 0. 


(23) 
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for all cotangent observables A, B. Because of this, we 
call (•, •)p the source metric, and (•, •)p, which can still be 
any of the contractive metrics, the target metric. This is 
the expression that we use to compute the effect of TZp. 
We can then compute the distinguishability via 

\\V\\lp = l3^{V,£TZpiV)fp. (24) 

For context, we note that if the both source and target 
metrics were the metric, then TZp would be a quantum 
channel, namely the transpose channel m, which is a 
useful approximate recovery channel in quantum error 
correction [HI dg. Classically, as both metrics reduce 
to the unique Fisher metric, TZp implements Bayesian 
inference from the conditional probabilities defined by £, 
relative to the prior p [g. 


III. GAUSSIAN STATES AND CHANNELS 

The theory of Gaussian states and Gaussian chan¬ 
nels can be formulated directly in an infinite-dimensional 
setting with uncountable number of degrees of freedom 
using the formalism of CCR (canonical commutation 
relations) C*-algebras. For simplicity, and clarity to¬ 
wards the intended audience, however, we do not work 
in full abstract generality. However, we use a formalism 
which should be relatively straightforward to generalise 
if needed. 

The beauty of the Gaussian formalism is that it pro¬ 
vides a one-to-one mapping between questions about a 
quantum system, to questions about a corresponding 
classical system. It is a mathematically rigorous formal¬ 
isation of the quantisation of free fields. 


A. CCR algebra 

Accordingly, we start by considering a classical phase- 
space, defined by a real vector space V equipped with 
a symplectic product (an anti-symmetric bilinear form) 
f^(/i ff)y /i 5 G One niay want to think of these phase- 
space points f,g as classical fields. The form cr(-, •) must 
be non-degenerate in the sense that if a{f,g) = 0 for 
all / then g = 0. For simplicity, it will be convenient 
to assume that V is equipped with a real scalar product 
{f,g) which is such that we can write cr{f,g) = (/, Ag) 
where A is an anti-symmetric real linear operator on V. 

Below, we will need to consider transformations of V 
corresponding to imaginary time evolution, which re¬ 
quires that we work in a complexification of V, which 
we call V^, where the scalar product is extended to a 
sesquilinear form via 

if + igj + ig) ■■= (/,/) + ig,g)-i{gJ) + i{f,g)- (25) 

The original V lives on as a subset of V'^, and we call 
/ S U a real vector of V^. 


A general classical observable is any real function on 
V, but we will focus only on observables which are linear 
functions on V, because the are the only ones which can 
be quantized unambiguously. An element / of the phase- 
space V can be mapped to a linear observable 4)/ defined 
by 4>f(g) := (f,Aq). We extend $ by linearity to the 
whole ofVc, . e„ := 4., + i#, if £ r 

The form cr(-, •), or equivalently A, defines the Poisson 
bracket on those linear classical observables via 

{$/,$g} = (/,A5)l, (26) 

for all f,g real, where !(/) = 1 for all /. 

The notation $j we use is there to evoke the fact that 
should be thought of as a smeared field observable. 
To make this clear, let us consider an example, where 
points of V are given by a pair of a field and its canonical 
conjugate: / = {<j), tt) G V. Here, the canonical field 
observable usually denoted “(('(a;)” would actually be the 
map on V which extract the component (^(x) G M from /, 
namely, the function / !->• (j){x). If the Dirac delta {6x, 0) 
was an element of V, we could write (j){x) = <l>(((5a;,0)), 
although typically it is not, which is why we need this 
somewhat more general formalism. 

Those linear observables are the ones that can 
be unambiguously “quantized”. More generally, though, 
one quantizes the bounded functions / i—> . Indeed, 

one defines the GGR algebra A{V) associated with V as 
that generated by the Weyl operators Wj defined by the 
relations WfWg = e~AfAg)y[r^_^^ and Wj = W-f for 
all f G V. One can show that the resulting C'*-algebra 
is essentially unique, and the operators Wf are unitary 
for all f GV. When extended to the complex V'^, these 
relations become 

WjWg = ei^^’^3^Wg_j, and wj = W_-j, (27) 

for all f,gGV^. 

If this algebra can be represented as that of bounded 
operators on a Hilbert space Ti, such that the unitary 
groups t !-)• W{tf) have generators, then these are un¬ 
bounded operators 4>/, such that 

Wf = e^^f, (28) 

satisfying the commutation relations 

[^f,%]=t(J,Ag)l. (29) 

Hence, if these operators exist, we can think of 
them as the quantizations of the classical observables 4)^. 
Moreover, like their classical counterparts, they are linear 
in their argument: 

^af+bg = ai>f + bi>g (30) 

for all f,g G and a,b G C. For most calculations, 
however, we only need to work with the Weyl operators 
Wf. 
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In the case of “second quantization”, such as non- 
relativistic quantum field theory, the classical phase space 
V is related to the Hilbert space H of “first quantized” 
wavefunctions as follows: the real and imaginary com¬ 
ponents of the wavefunction ip play the role of canonical 
conjugate variables. Hence V is just H conceived as a 
real vector space. If (V’l^) denotes the complex scalar 
product of "H, we may use on V the real scalar product 
('i/', (p) := Re {'iplp). The symplectic form is given by the 
linear operator A = —il as 

^(■0, «<>) = ('*/'. A0) = Im(V'|(/>). (31) 

In this example, the complexification V'^ is not equal to 
TL (it has double the dimension). Consequently, we avoid 
this formalism. 


B. Gaussian states 


It can be deduced from Equ. (27) that the whole alge¬ 


bra of observables A{V) is linearly spanned by the Weyl 
operators Wf. Hence, a state p is entirely characterised 
by its value on those. These values are summarized by 
the state’s characteristic function f i— p(Wf). 

A Gaussian state is one whose characteristic function 
is Gaussian: p{Wf) = e~^AAf)+i-(fJo) for all / G H and 
some /o G V, where A is a symmetric bilinear operator 
on V. When extended by linearity on V'^, we therefore 
have 


= A = A. (32) 

Example of Gaussian states are the thermal states of 
free bosonic held theories. A similar formalism exists 
for fermionic helds, but we treat only the bosonic case 
in this paper. In what follows, we assume /o = 0 for 
simplicity, hence 


p{Wf) = e 2 


_ p-|(/A/) 


(33) 


for all / G V'^. This equation together with Equ. (27) 


allows one to compute the expectation value of any op¬ 
erator. Eor instance, the expectation values of products 
of smeared held operators can be evaluated by successive 
differentiation of this expression. We hnd 


= (/,(A+lA)g), 


Hence the operator A determines (and is determined by) 
the real part of the 2-point correlations functions. 

The above expression also implies that the complex 
operator A+ must be positive. 


The operator A + may be thought of as the Hamil¬ 
tonian (rather than Lagrangian) version of the thermal 
propagator. Here, we will call A the covariance operator 
as is traditional in the Gaussian formalism. 


1. Classical Gaussian states 

The above formalism is almost identical for classical 
statistical theories. The only difference is that A = 0, so 
that the held operators are commuting. Indeed, the 
Poisson brackets plays no direct role at this level. The 
Gaussian states can then be interpreted as thermal states 
of free classical helds living in the phase space V, for 
some Hamiltonian related to A. Specihcally, given the 
quadratic classical Hamiltonian / i—)■ i(/, 77/), we hnd 
that the corresponding thermal state at inverse tempera¬ 
ture /3 is the Gaussian dehned by the covarience operator 

(36) 


C. Gaussian channels 

A channel in quantum theory refers to the most gen¬ 
eral map from states to states which is consistent with the 
probabilistic interpretation of the convex combination of 
states, as well as with the system being part of a larger 
one. In the Heisenberg pictures, it must be a linear map 
between algebras of observables that is completely posi¬ 
tive, and which preserves the identity. Glassically, these 
map correspond to all stochastic maps. In quantum the¬ 
ory, channels include unitary transformation, but also 
maps which add noise to the system, and correspond to 
the evolution of open quantum systems. We will not go 
here into details of this definition, because we will simply 
focus on a special class of channels which map Gaussian 
states to Gaussian states, the so-called Gaussian chan¬ 
nels. These also include all unitary transformations as a 
particular case. 

As with states, channels are characterised entirely by 
their action on the Weyl operators. We consider channels 
S with the following action 

£{Wf) = Wx/e-5(7A/) (37) 

where X and Y are linear on V, i.e., real on V'^\ = 

and Y'^ = yl. We can easily check that this maps 
Gaussian states to Gaussian states, provided certain con¬ 
ditions on the linear operators X and Y. li p a is the 
Gaussian state defined by the covariance matrix A, we 
obtain using Equ. ( [33| that 

PA{S{Wf))=pxtAX+Y{Wf). (38) 

Erom the requirement that X"^AX -\-Y -|- |A > 0 for all 
A -I- IA > 0, we obtain the condition 


A+'-A>0, 


(35) 


Y - -X^AX + -A > 0. 
2 2 ~ 


(39) 
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We note that classically, with A = 0, this simply re¬ 
duces to the condition F > 0. 


IV. DISTINGUISHABILITY NEAR GAUSSIAN 
STATES 

We now want to obtain the adjoint map TZp using 
Equ. ( [^ when p is a Gaussian. The strategy is to eval¬ 
uate all components in terms of the Weyl operators W/, 
or more conveniently, in terms of the functional 

Gf (40) 

Since those operators are not self-adjoint, we need to ex¬ 
tend the metric to all operators, making it sesquilinear. 
Using Equ. 0 , this is 

{A,B)l:=p{A^ep{B))- ( 41 ) 


Moreover, the fact that it is canonical means that it pre¬ 
serves the symplectic form: {Rf)"^ARf = A. Another 
important property of these operators is R^ = from 
which we obtain that 

iA+^-A)Rf = iRf)\A+^-A). (49) 

Also, by computing the components of the map G/ i—>■ 
pGf = Gj^aj^p from both expression, we obtain 

(A+^-A)Rt = A-^-A. (50) 

Explicitly, the Kubo-Mori metric is 

{Gf,Gffp=[ ( 51 ) 

Jo 

Then the definition of TZp in Equ. ([2^ becomes 


We will make use of the fact that, assuming A is the 
covariance matrix of p, 


p((G^)lG^) 


(42) 


for all f,g G U"'. This can be computed directly using 
Equ. (27) and Equ. (33). 

Abbreviating 


B = X^AX + Y 


(43) 


for the covariance matrix of the coarse-grained state po£ 
(See Equ. ([38|)), we also have 


- ^(/.(B+tA)g) 


pi£{{GfyG^))=e 


(44) 


In terms of those functionals, Equ. (38) becomes simply 

£{Gf) = Gif. (45) 

The particular metric £lp defined in Equ. ( |17[ ), corre¬ 
sponds to the operator monotone function 


6^{x) = 


X —1 
logx 




(46) 


Recall that 0p is given by applying 9 (through func¬ 
tional calculus) to the superoperators of imaginary time 
evolution, namely the transformations X 1 —>■ pXp~^ = 
Xq 0 h_ For a Gaussian state this generates a group 
of (complex) canonical transformations which can be rep¬ 
resented by linear operators Rf (where A stands for the 
covariance matrix defining p) on the complexified phase 
space: 

= WjiAf. (47) 

Since p is invariant under the imaginary time evolution 
it defines, we have R^ARs = A, hence 


(7^p(G^),Gf)po£ = {Gf,£{Gf)fp = {Gf,Gip); 


= r ^{RtfAA+iA)Xg)^^^ 
Jo 


Equ. (52) can be used to prove the following. Let 
be the linear space of polynomials obtained by differen¬ 
tiating G^{f) with respect to its argument n times, i.e.. 


n := 






dti . . . dtr, 


■ fl, ■ ■ ■ , fn G V 


} (53) 


where the subscript 0 indicates the the derivative is eval¬ 
uated at ti = • • • = tn = 0. These are polynomials of 
order n in the field operators. Moreover, the span of 
to is the space of all polynomials in the fields of order 
n (or lower). 

Because the exponent in the last term in Equ. ( |5^ is 
linear in both / and g, the whole expression is equal to 
zero whenever one differentiates it with respect to / and 
g a different number of times (at f = g = 0). Moreover, 
Equ. (51) shows with the same argument that is or¬ 
thogonal to V'y in terms of the Kubo-Mori metric at p 
whenever n y m. 

Similarly, if we assume that the target metric is defined 
by an operator monotone function of the form 9{x) = 
f x‘^dp{s), for some measure p, which includes the Kubo- 
Mori metric, Bures metric and metric, then we have 
explicitly, 

{Gf,Gf)po£= j (54) 


which implies that V'^ and are orthogonal also in 
terms of that target metric. Taken together with the 
completeness of all these polynomials, this proves that 


^-spHGj^spH ^ Q 


A 

Rff 


(48) 


np{vi) C VI- 


(55) 
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Since also, for Equ. (451, £{1^%) C "P^j we obtain that 
for any Gaussian channel £, and any Gaussian state p 
with covariance operator A, 


£n,{vi) C VI. (56) 


Moreover, since is orthogonal to whenever n ^ m 
in terms of the Kubo-Mori metric (•, •)p, then 

(K,£7^p(Kr^))^ = 0 (57) 

for any 14, G with n ^ m. This implies that we can 
restrict the problem of computing the components of the 
linear map TZp to each subspace V^ independently. 

With some extra assumption on p and the channel, the 
same argument can yield a more detailed result which 
will be useful below. Suppose that X, Y, A and A are 
all jointly block diagonal for some decomposition V = 
14 of the classical phase space (assuming that k is 
an integer for convenience). Let us define the operator 
spaces 


V. 


ki<‘“<kj 






dti . . . dtr, 


■ h&Vk^i = 1 


... ,n|. 


(58) 


Then, following the same argument as above, we see that 
these spaces are all orthogonal to each other in the Kubo- 
Mori metric at p, and also that 


(59) 

An example is the Klein-Gordon example studied be¬ 
low, where k index momentum modes, and the spaces 
'^ki<...<k finite-dimensional, allowing for an exact 
solution for each family of modes. However, for some 
specific metrics, or with a rather innocuous simplifica¬ 
tion, the solution can be made much more explicit. 


A. Exact solution for metric 


It follows that 


STZpiGf) = f £{Glf)ds = Gij^jds, (64) 
Jo Jo 


and 


' p 


{Gj,£np{Gf)fp= [\s{Gj,GijJ; 

Jo 

= (65) 

Jo 




where 


P={A+ + ^A)-iAt(A + ^A). (66) 

Note that = P because 

R%^{B + ^A)-i = R\/^{B + ^A)-i(i?f,/jt, (67) 

One can then compute ||K|||_p for any operator V by 
differentiation of the last expression in Equ. (651 with 
respect to / and g. In general, one can write the operator 
V as 


V = DjGf (68) 

where Dj is a differential operator with respect to the 
observable / (we give some explicit examples below). Us¬ 
ing Equ. ([T^, this yields 



where 


The quantity IjUlj^^p, or more generally all components 
of the coarse-grained metric, can be computed exactly in 
full generality for the so-called metric, defined by the 
function 

0{x) = , (60) 


which yields 


(Gf,Gf)pof 


(61) 


Indeed, we see from comparing this to Equ. (52) that 
TZniGf) = fn G?,.ds provided that h{s) is such that for 

all0<4<4„dfor.ll/,Ser, 


(Rf/Ms), (B + JA)g) = (Rff, (A + ^A)^^), (62) 

namely h{s) = Jsf where 

J. := R^,/^iB + '-A)-^XHA + ^-A)Rt (63) 


P{s,t) = iRf)^PRt (70) 

The result from the action of these differential operators 
on the exponential generator can be expressed in terms 
of familiar Feynman diagrams, with the propagator given 
by relevant components of P{s,t). 

For instance, a differentiation of order n yields a Feyn¬ 
man diagram with n edges in total, which represents a 
contraction of the operator 

1 

P{s,t) ^ ^ Pis,t) ds dt. (71) 


B. Exact solution for the classical Eisher metric 

All quantum contractive metrics reduce to the Fisher 
information metric when all operators commute. For our 
purpose, this implies A = 0, as well as 7?^ = 1 since the 
imaginary time evolution is trivial: p~'^Xp‘^ = p~''^p^X = 







X for all X. Hence we can use the above result to directly 
get the classical solution simply by setting A = 0 and 
i?i = 1 in Equ. (661, yielding 


P = AXB-^X^A = AXiX'^AX + vy^X'^A 
= A{A+{X'^)-^YX-^)-^A, 

and 

(G/,£7^p(G^))p = e(^’^9), (73) 

where we removed the S label since all metrics used are 
the same Fisher metric in this case. This result is equiv¬ 
alent to that obtained in Ref. [4]. 


C. Approximation for large noise 

There are contexts where the channel, defined by X 
and Y, is noisy enough that the coarse-grained propaga¬ 
tor is approximately independent of the imaginary time 
s: 


(H + cs (74) 

for any s € [0,1]. We will see below an example where 
this is the case. In this case, no matter what target metric 
we use, 

(75) 

so that we can use the same method as in Section flV Al 
to obtain 

(Gf,£7^p(G^))® ~ f (76) 

with 

P = {A + ^A)Aif-pipXt(A + ^A). (77) 


V. EXAMPLE: KLEIN-GORDON FIELD 

The classical Klein-Gordon Hamiltonian can be writ- 


ten as 



1 

^=2 

J dx {nixf + (l){x){ni^ - d'^)^{x)), 

(78) 

where tt(x) and (pix) are canonical conjugates, 
canonical change of variable 

The 


f dx[cosikx)^{x) - sin(fcx)^(x)] 

(79) 

n, = j 

f dx[oJk sm{kx)^{x) + cos(fcx)ft(a:)]. 

(80) 

(81) 


yields the decoupled form 

77 = - y dk{tv\ + . ( 82 ) 

where = yjk"^ + rrP and the Poisson bracket is 

= S{k - k')l. (83) 

Instead of attempting to quantize this directly, we con¬ 
sider the discretisation 

^ = 2 ^ “ 2 (84) 

k k 

where 

$fc = L“5$fc and H^ = L'^Hfe (85) 

satisfy 

{$fe,nfe} = i. ( 86 ) 

The “infrared regulator” L‘^, where d is the dimension of 

space, can be thought of as a volume. The original theory 
with continuous momenta k is recovered for L —^ oo. 

The discrete Hamiltonian simply represents a discrete 
set of decoupled harmonic oscillator. Therefore, the cor¬ 
responding quantum Hamiltonian is 

H = '^Ukalak (87) 

k 

where the sum is over some discrete set of modes fc, and 
[afe,aj,] = 1 and [ 0 ^, 0 ^/] = 0 for k ^ k'. The quantum 
versions of the observables (j)k and tt^, which we denote 
by the same symbols, are 

= ^yl/2oJ{al + Ofe), 

Hfc = i^/^JJ2{al - Qk). 

Since there is no interaction between modes fc, the 
Gibbs state is of the form p ~ 0^, pk- In the basis com¬ 
posed of the classical observables ($fc,nfc), the symplec- 
tic form is 

Afe = . (88) 

Using the creation and annihilation operators, one can 
easily find the expectation values of products of two field 
operators, and hence the component of the covariance 
matrix A^ defining the state pk'- 

A, = lcoth(^)(^-Q'c (89) 

One can obtain the components of Rf‘‘ by solving the 
imaginary time equations for the harmonic oscillator, 
which yields 

nAfc _ ( cosh(^a;fcs) -iujk sinh(/3a;fcs)\ , . 

® yismh{f3ujks)/ujk cosh(/?Wfes) J ' ^ 
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We could proceed using the real phase space coordi¬ 
nates and Ilfc. However, the coarse-graining channel 
that we will use takes a simpler form in terms of the 
complex variables 


But recall that X and Y must satisfy Equ. (39) in order 
for po£ to be a valid state, or, said differently, for £{A) 
to be a positive effect whenever A is. On mode k, the 
equation reduces to 


4’k = 2 ~ 

n_fc) — — ‘b-fc)- (92) 

(93) 

These observables are the standard Fourier modes used 
in scalar field theory. They are related to the original 
fields and 7r(x) simply through 


W + (l-e-'=''^')^A>0. (101) 

For k l/cr, this just implies that Yk > 0. But for 
k ^ 1/(7, this yields the non-trivial relation 

ylvl > 1 - ( 102 ) 

This can be understood intuitively by looking at the co- 
variance matrix 


(/(x) 

= hm L5 

L—¥00 

(94) 


k 


7r(x) 

= lim P^^-Kk- 

L—¥oo 

(95) 


k 

(96) 


Together with (^_fc,7r_fc) = this is just a 

complex change of coordinate on the four-dimensional 
subspace of corresponding to modes k and —k. Re¬ 
call that A and A both define sesquilinear forms on V'^. 
Therefore, if we denote the components of this coordinate 
change by the four-by-four matrix P, Ak © A-k trans¬ 
forms to P^Ak © A-k)P in the new coordinate system. 
One can check by direct calculation that they have the 
exact same form as before: P^Ak © A-k)P = Ak® A-k 
and Pl(Afe © A_fc)P = A^ © A_fe. 

From now one, by “mode k” we mean either the sub¬ 
space of phase space spanned by {4>k,T^k), or the corre¬ 
sponding subsystem in the quantum theory. 

In order to proceed further, we need to fix a channel £ 
via the operators X and Y. We use the linear operator 
X defined by 


X^k = (97) 

XiTk = (98) 

(99) 

The operator Y defines a sesquilinear form, which we 
take to be block-diagonal Y = ^keM where Yk can 
be represented by the matrix 



This is a variation of the channel used in Ref. [J] , but this 
operator Y is different as it couples the real modes k and 
—k when expressed in terms of the coordinates [4>k,T^k)- 
The parameter a characterises the maximal precision at 
which space is resolved by the observer, and the values 

d d 

L'iy^ and characterise the precision at which the 

field and canonical field values (/(a:) and tt(x) respectively 
are resolved at each point in space. 


^k — x\AkXk + Yk 


(103) 


which defines the state £k{pk)- When Xk — 0, the covari¬ 
ance matrix is just Yk^ hence it must satisfy the Heisen¬ 
berg uncertainty relation given by Equ. (102). 


It is possible to write Bk in the same form as Ak, but 
in terms of redefined frequencies and temperatures: 


Bk = i coth(%i) {< ® ) . (104) 

\ U ^kJ 


Let us define 


vl ■■= yl+^k\ coth ('^''2 ")e ^ . 

(105) 

(106) 

Then we have 


, Vk 

UJu — — 

Uk 

(107) 

and 


/?(, = —2coth“^(2ufcUfe). 

Vk 

(108) 

Moreover, assuming y^y^r ^ 1, which implies 
for all k, we obtain 

UkVk » 1 


P'k^ 


(109) 


With these definitions, Rf'‘ is given simply by substitut¬ 
ing LOk —t w/ and ,0 —)■ /3(. in Equ. (901. 


But since then ~ 


<C 1, each mode k is 


in the high temperature limit where all quantum metrics 
reduce to the classical one. Indeed, we obtain simply 


Rf ^1 + Oiy^y^)-^ (110) 


for all s. This puts us in the situation described in Sec¬ 
tion |IV21 with 


l^app — 

k 


{1 



( 111 ) 
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As an example let’s compute the distinguishability of 
(j)k and TTfe, which are given by the first order derivatives 
of G^. Differentiating Equ. (76) once with respect to / 
and g, we obtain that 


Mfc := 


/'{(j)k,STZp{(j)k))p {(j)k,£'R-piTrk))p 

\{nk,£'Jlp{cj)k))p (7^fc,£7^p(7^fe))^ 


= Kl 


where 




= f\Ak + 

Jo 


0 \ -k^a^ 

-2 ) Kse 




( 112 ) 


0 


(113) 


Hence we obtain from the diagonal components of Mk 
that, with ^ 1, 


UkWlp^ 






4,,2 


and 


ujlu- 


fc“fc 


\\T^k\\s,p - 


(114) 


For the calculations involved higher order polynomials, 
we avoid doing the integrals over t and s by considering 
only the zero temperature limit j3 —>■ oo. Recall that, at 
finite temperature. 


D«(F) = -/3 


i d e 


-P(H+tV) 


dt 


|t = 05 


(115) 


for V self-adjoint and of zero expectation value. At zero 
temperature, the state /Z\t=Q should become 

the projector on the ground state \£lt) of the quadratic 
Hamiltonian H perturbed by tV. Hence, if we write 
|0() = |Do) -I- t|Di) -I- we have 


hm /?(lE,l/)^ = -^Tr(lE|D,)(D,|)|,=o 
P^ao ^ dt 


= -(Do|lE|Di)- (Di|lE|Do)- 

Moreover, from perturbation theory, 

\ni) = -H-^v\no), 


(116) 


(117) 


where H is the Hamiltonian, shifted so that the ground 
state energy is zero, and H~^ is defined to be zero on the 
ground state. Hence, we can rewrite the Kubo-Mori met¬ 
ric at zero temperature, now for the complexified version, 
as 


li^J{Gf,Gffp 

= p{{Gf)^H-^Gf) + p{GfH-\Gf)^). 
For our example, using the Hamiltonian, 


H 


E 


0 Jka\ak 


(118) 


(119) 


we can obtain the components of the metric in the sector 
of distinct modes fci,..., at /3 = oo. From differentiat¬ 


ing Equ. (118), using Equ. (42) and Equ. (51), we obtain 


that 

^Ihn /3 I ^{Ak^+^A)R, 


nl n 


i=i 

1 




(g)(A^ + -A)-h(g)(A^--A) 

i=i i=i 


( 120 ) 


where Af = lim^^oo Ak = \\‘^’^i ^ . 

VO 

We used the fact that 
dt,--- dt^ds, - - - 

= '^ifki, Pgk ^,) • • • (/fc„, Pgk„^), 


( 121 ) 


where the sum is over all permutations tt of {1 ,..., n}. 
From Equ. (76) and using the above results, we can 


now obtain any component of the coarse-grained met¬ 
ric. For instance, in distinct modes fci,..., kn, if we de¬ 
note the components of the Kubo-Mori metric (given in 
Equ. (120)) by then the components of the 

coarse-grained metric in our approximation are 


^ki,...,k„ 


i=i 


( 122 ) 


We consider first the example of the mass term in the 
Hamiltonian, namely 


iV = ^ / (j){xYdx — cl 


(123) 


where c G K is there so that {V 2 )p = 0. In terms of the 
variables (pk, this is V 2 = \ '^f.(j)k4'-k — c'l. In terms of 
the generating operator Gj , this is just 




a" 




2 ^dsdt 

k 


t=s=0’ 


(124) 


Indeed, the identity component is automatically absent 
from this expression due to the fact that polynomials 
of different orders generated by G^ are automatically 
orthogonal in terms of the Kubo-Mori metric, and the 
fact that (K, l)p = {V)p. 

Let’s use the shorthand P* := where P is 

given by Equ. By applying the above differentiation 
to both / and g in Equ. (76), and using Equ. (121), and 
the fact that {(j)k,Ps4>k') = 0 ii k ^ k', we obtain (still 
for 2 / 0 ?/,r » 1), 


I Eg 


k ' 

Ei: 


ds dt {(j)k, PgVfe)(</>-fc, Pl(l}-k) 

ds dt[(j)k ® (j)-k, {Ps ® Ps)(l^k ® (t>-k)- 

(125) 
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Using Equ. (120), each term of the sum can be com¬ 


pute explicitly in the zero temperature limit. The distin- 
guishability density of V 2 , which we denote 


d(V 2 ) = lim L-‘^\\V 2 
L—^oo 


\Ip 


(126) 


is then obtain simply by replacing the sum by an integral 
in Equ. (125), which is 


d{V2) ~ 


1 


dk 






1 

~~d 4 


1 


-'fc j 


(127) 


where Uk and Vk are given by Equ. ( |105[ ) and Equ. (106 1 , 
but in the limit /3 —)■ 00 , i.e., with coth(^^) —)• 1 . 

We also want to consider the quartic interaction term 


^^ = 4! 


11 m‘ 


dx - 


(j){xYdx — cl, (128) 


where the identity component compensates the opera¬ 
tors’ non-zero expectation value, and the second term is 
a counter term obtained from renormalisation (only to 
first order in perturbation theory). It is needed because 
the quartic term alone yields expectation values which 
diverge in the continuum limit (when an ultraviolet cut¬ 
off is removed), and hence does not by itself constitute a 
valid perturbation of the Hamiltonian in the continuum. 

An ultraviolet regularisation is given simply by 

This is due to the orthogonality of the polynomials gen¬ 
erated by in terms of the Kubo-Mori metric, which 

guarantees that adding V 4 to the Hamiltonian has no 
influence on moments lower than four (to first order in 
the coupling constant). In particular, this means that V 4 
does not influence the value of the “macroscopic” mass 
that is defined from the second moment in (/)(x), which 
is the role played by /r in Equ. (128). 


To proceed, we need to express this interaction in 
term of our variables 4>k- Observe that / (j){x)‘^dx = 
f dki... dki ki) (j)k^ ... We use the infrared- 
regulated operator 


1 


^4 = ^ E 


. ,^4 


dti ■ ■ ■ dti 


(130) 


which yields at zero temperature 

1 


d{V4) ~ 


4! 26 J 

-cr^Sikf r 


Ln, 


dki ■ ■ ■ dk^ 5{Tiiki) 

1 1 


i '^ki 


U.niojl 


E 


U ^ki 

1 


(132) 


222222 
uf uf vf vf LOT ojf 

fevrCl) ^^(2) fc7r(3) fc7r(4) fejr(l) ^^^ 2 ) 


where the products are over i = 1,... ,4 and the last sum 
is over all six permutations tt of the set {1, 2, 3,4} which 
are such that 7r(l) < 7r(2), and, of course, Uk and Vk are 
taken at B = oo. 


A. Classical version 

For comparison, we also consider the same calculation 
but for the classical Klein-Gordon field, that is. 


o' 


A, — 

0 


(133) 


The operator P from Equ. (72) is P = ®kPk with 

Pk = Al{Ak + Yke^"’^")-\ 

In particular, 

1 


(134) 


{(pkj Ppk) o 2 I 02 4 2 k^c 

Bui + 


We obtain 


and 


d{V2) = - dk{Bk,PBk)\ 


(135) 


(136) 


17 4 

diV/B) =j dki" ■dki5{Y^ki)'^{4)ki-,P(t>ki)- (137) 

This is the value of the “basketball” Feynman diagram, 
albeit for a modified propagator. A partially analyti¬ 
cal solution for the standard propagator can be found 
in Ref. [17) . However, one can use a more versatile nu¬ 
merical method to evaluate it, such as explained in Sec¬ 
tion |VE 1 


so that 14 = limi_,.oo P‘^14- Applying these derivatives to 


Equ. (76) like in the previous example, we obtain (still 
for » 1 ), 


17 ^ 

II^4|IL=^4! E dsdtdiY,kBll{Pk.,PlPkB 

k"..k.d i=i 


( 131 ) 


B. Comparison with regulated Bures metric 

Other authors have considered the raw Bures metric 
dependent on a sharp momentum cutoff, e.g., in Ref. [2 
[2]. This is equivalent to using a coarse-graining channel 
S defined by E = 0 and X the projector onto modes |fc| < 
1/cr for some scale a. Technically, the target algebra of 
£ should not contain any mode larger than 1 /cr, or this 
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would violate Equ. |39[ hence X is actually an isometry 
(satisfying X^X = 1). The resulting channel £ then 
simply performs a partial trace over all modes |A:| > l/tr. 

In this section, we want to examine under what con¬ 
ditions such a quantity would match, possibly approxi¬ 
mately, the coarse-grained metrics we computed. 

The “raw” distinguishability of a Hamiltonian pertur¬ 
bation in the Bures metric {X,Y)p = Tr(XH“^(y)) is 
dehned by 

^p{^) = + P^)- (138) 

In th e cas e where p is pure: p = |Ho)(Ho|, we can use 
Equ. ( |1I8[ ) to relate it to the Kubo-Mori metric as 

\\mJ{Gj,G^)l = 2{H-^Gj + GjH-\G^)p. (139) 


What we want to compute is simply our coarse-grained 
distinguishability when £ = id. Note that this does not 
imply TZp = id because the “source” metric (Kubo-Mori) 
is different from the target one (Bures). Instead, we have 

^Ihn /3(7^(G^),G^), = \\^J{Gj,GX 

= 2{H-^Gj + GjH-\Gj)p. 

(140) 


Hence, 

n{Gj) = ^{GfH-^ + H-^Gf) + G((l//3)°). (141) 

Using also Equ. ( |118[ ), the distinguishability is then 
{Gf,£n{Gf)fp=^m2l3{Gf,GfH-^+H-^Gffp 
= 2p{{GfyH-^Gf) + 2p{GfH-^{Gfy). 

(142) 


Hence components of this metric in distinct modes 
ki,... ,kn are given by the matrix 


2 


i=l 



(143) 


In order to compare to the components of the coarse¬ 
grained metric in Equ. (122), observe that, for the Klein- 
Gordon field state at /? = oo. 


(Hfe ±-A)Hj, =F 2^) = 0 (144) 

{Ak ± ± ^A) = 2{Ak ± '-A). (145) 


VI. APPLICATIONS 
A. Wilsonian relevance 


An interacting QFT comes with a regulator such as 
a momentum cutoff. In general, even if it may loosely 
correspond to the energy beyond which the theory loses 
validity, it is not considered to be physical, or at least its 
value is not known. As such, the predictions of the theory 
at low energies must be independent of it. In general, 
the parameters of the theory (entering the Hamiltonian) 
must be made dependent on the cutoff so as to yield 
constant predictions. This function of the Hamiltonian 
on the cutoff is the renormalisation group (RG). 

Wilson argued that the RG can also be reinterpreted 
in a different way [18]. As the energy cutoff is lowered, 
the theory may simplify in the sense that some “coupling 
constant” in the Hamiltonian, such as the parameter A in 
H + AU, tend to zero. Wilson interprets this as meaning 
that the term V is not needed to make correct predictions 
at low energies. Hence the simpler theory defined by 
the Hamiltonian H is a good effective description of the 
system at low energies. 

A more detailed way of addressing the question of 
whether V is detectable under certain experimental con¬ 
straints is of course the calculation of an actual statistical 
distinguishability between H and H + XV as a function 
of explicit resolution parameters. 


Specifically, the quantity d{V) defined in Equ. (126), 
tells us how easy it is to distinguish a Gibbs state for 
H compared to one for H + eV, to lowest order in e, as 
a function of the various resolutions cr, j/^, per unit 
of volume. The behaviour of d{V) as a function of the 
“physical” resolution parameter a should tell us directly 
how relevant U is: if it decreases as a increases then one 
may deem it “irrelevant”. But this doesn’t quite work, 
because, since all distinguishability measures are contrac¬ 
tive, d{V) always decreases as a increases. However, it is 
a density: it corresponds to the effective distinguishabil¬ 
ity of the term V for an observer having access only to 
some fine volume. When the resolution cr gets worse, the 
observer has effectively access to fewer “pixels” per unit 
of volume, which is why d(V) decreases. 

To meet the standard concept of relevance, we should 
rather ask how the distinguishability of V depends on a if 
the observer still has access to the same volume relative to 
a. This means considering instead the unitless quantity 


^(I/) := a^d{V), 


(146) 


Hence, the raw Bures metric is given by the substitu¬ 
tion — )■ vl —>■ ^ujk and X —>■ I in Equ. (122). 

This suggests that, for the Klein-Gordon example, our 
coarse-grained metric essentially matches the raw met¬ 
ric (with a momentum cutoff) at zero temperature when 
uj. ^ and v^. Wk- This would violate the un¬ 
certainty condition 1, but this condition is not 

fundamental as it can be easily alleviated by making 
and j/tt dependent on k. 


which can now increase or decrease. This extra factor 
a‘^ serves the same purpose as the active rescaling of the 
fields in Wilson’s approach to renormalisation. 

More specihcally, we are interested in the polynomial 
behaviour of 5{V) as a function of a. Hence we can di¬ 
rectly extract the degree of the polynomial as 


a(U) := 


91og5(U) 
9 logo- 


(147) 
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For instance, consider the case where yrj, is small, which 
also requires large given the uncertainty relations. We 
have, roughly, 

d(F2)oc [ dk\ (148) 

J\k\<l/a 

Therefore, for cr ^ ^ (or a massless theory), d{V 2 ) = 
0 (0-4-^), and for ct > ^, d{V 2 ) = ©(ct"^). This yields 
a{V 2 ) — 4 below the mass scale and a{V 2 ) — 0 above the 
mass scale. 

This can be understood as follows: as the “pixel size” 
cr increases, more information about the mass is gained 
per pixel, until they reach the mass scale (which is also 
the correlation length), at which point no more new in¬ 
formation is gained. Hence Information about the mass 
is found in large scale features of the state. 

This corresponds to the Wilsonian analysis which 
would say that V 2 is relevant for a massless theory, and 
marginal in a massive theory. 

With a large field uncertainty the scaling is quali¬ 
tatively different. Indeed, below the mass scale, 

d{V 2 )^ [ dk^, (149) 


a 



FIG. 1. Scaling exponent a for the distinguishability of the 
(j)^ interaction as a function of the scale (resolution) a, where 
CTO := 1/m is the correlation length. For both curves, = 
10 ^°/cto, and the dimension of space is d = 4. The solid curve 
is for y'^ = 10“®cto (a valued marked on the scale axis by 
the vertical dotted line). For comparison, the dashed curve 
is for the much smaller = 10“^cto, for which a behaves 
essentially as predicted by Equ. ( 1501, i.e., a ~ —2 below the 
mass scale. 


The last term involves an expectation value in terms of 
the probability density K 1 — f{K)/ f dKf{K). Defining 

H{K) =-log f{K), (153) 


which yields a{V 2 ) — 6. 

As a further illustration, we want to compute the dis¬ 
tinguishability scaling a{Vi) for the 4 >‘^ interaction term 
from Equ. (132). We can make a similar “back of the 
envelop” analysis for the behaviour of ^(14) for y^^ very 
large and y^ vanishingly small. This is 


At TO = 0, Wfc = k. Counting the powers of k inside the 
integral yields the guess a(14) ~ 2(3 — d) for 2 < d < 6, 
which can be verified using the numerical method ex¬ 
plained below. Hence we recover the expected marginal 
dimension dp = 3 above which the interaction is irrel¬ 
evant at large scales. However, various other results can 
be obtained depending on the field resolutions y^j) and y,r, 
as can be a nticip ated by the resulting different powers in 
ujk in Equ. 

In order to obtain an efficient numerical method to 
evaluate the integrals, consider that, if 


d(H) = j dKf{K), (151) 


where K = (ki ,..., kn), then 
i91ogd(y) 




cHogcr 
1 


^ / dlogf{K) \ 

\ da If 


(152) 


the integral can be evaluated using the metropolis algo¬ 
rithm for the classical Hamiltonian K 1 —>■ H{K) at tem¬ 
perature 1. 

If we keep y^ and independent of ct, those values 
introduce two lengthscales in additions to the mass, at 
which the behaviour of d(V 4 ) changes. The effect of an 
increase in y^ is shown in Figure]^ In this example with 
four spatial dimensions (d = 4), an increased value of 
y^ renders the interaction V 4 relevant {a > 0) instead of 
irrelevant (a < 0) below the mass scale for a < y^. 


B. Dimensionality reduction 

An approach to renormalisation proposed in Ref. HIS] 
consists in determining the eigenvalues and eigenvectors 
for the linear map STZp, which is self-adjoint with re¬ 
spect to the source metric. The original interpretation 
of the calculation was different, because TZp (which was 
then called TZj^) occurred geometrically as the pull-back 
induced by £ on cotangent spaces to states (we refer to 
Ref. H for details). However, technically the calculations 
are identical. 

In that approach, the eigenvectors are interpreted as 
observables, and the eigenvalues (which are always in the 
interval [0,1]) tell us how much these observables lose 
“distinguishability” under coarse-graining. A smaller ef¬ 
fective tangent space (to states) is given by ignoring those 
observables which lose too much distinguishability, i.e., 
which are not “relevant” enough, by deeming two tan¬ 
gent vectors effectively identical if they give the same 
expectations values for the more relevant observables. 
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For the Klein-Gordon field example, the components 
of £TZp as a linear map can be determined using the ob¬ 
servations made at the end of Section EYi namely that 
STZp is block-diagonal, with each block associated with 
momentum modes ki,... ,kn- For yT^y^ S> 1 and at zero 


temperature, one can directly see from Equ. (122) that 


the components of one block are given by the matrix 


'u7^ 0 


1=1 


0 V 


-2 


K 


k^cr^ 


(154) 


where the matrices are equal to the right hand 

side of Equ. (|120|). 


For instance, for n = 1 one recovers the single-mode 
result presented in Ref. [1] (at /3 —>■ oo). For n = 2, and 
ki = —k 2 = k for simplicity, we obtain the eigenvalue 


m = 


g-2fc^o'^ 

2PuJkulvl 




with corresponding eigenvector (observable) 
Al = (j)kT^-k + T^-k4>k 

and 




= 


4/3 




W/c 


o(r"). 


for 


A2=uiu:l\^k?-vt\cl>k?. 

The other two eigenvalues are of order /3“^. Two orthog¬ 
onal polynomials spanning that space are given by 

^3 = *(</'fc7r-fc - 7r_fe(/)fc) and A^ = ujl\(j)k\^ + \'Kk\^■ 


VII. OUTLOOK 


We developed a practical way of computing the ef¬ 
fective distinguishability between Hamiltonian perturba¬ 


tions, taking experimental limitations explicitly into ac¬ 
counts. These resolutions parameters are needed to ob¬ 
tain a finite measure of distinguishability. We developed 
a calculation procedure applicable directly to Gaussian 
states (free fields), to arbitrary perturbations, and for all 
contractive metrics. Moreover, we showed in the context 
of the Klein-Gordon field that all contractive metrics give 
the same result provided that the phase-space resolution 
is much coarser than h. 

Although the quantities computed, namely the compo¬ 
nent of a metric, only give the distinguishability between 
a state and a perturbation of the state to first order in 
the coupling constant, this can be used in principle to 
compute geodesic distances between any two state (or at 
least upper bounds to it). This is where the geometric 
nature of the formalism becomes useful. For instance, 
the distinguishability between a state pab on two sys¬ 
tems A and B and the product state pa® Pb (where 
PA and pb are reduced states of pab) is a genuine mea¬ 
sure of entanglement when pab is pure, such as a ground 
state. This can be used to compute an upper bound 
to the ground state entanglement on any region A. For 
instance, if p is the ground state of a local field Hamil¬ 
tonian, then a path to the product state pA ® Pa'^ may 
be obtained by progressively increasing the field’s mass 
on the boundary. In the purely Gaussian case, the re¬ 
cent results in Ref. m may provide a more direct route. 
However, our method can in principle be extended per- 
turbatively beyond Gaussian states using the approach 
proposed in Ref. |3] in the classical case. Namely, a per¬ 
turbation of a Gaussian state p yields a perturbation of 
the linear map STZp, whose components can in principle 
be computed to each order. 

A generalisation to fermionic fields should also be pos¬ 
sible using the fermionic Gaussian formalism, such as de¬ 
veloped in Ref. [20] . 
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